Specific heat across the superconducting dome in the cuprates 
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The specific heat of the superconducting cuprates is calculated over the entire phase diagram. A 
d-wave BCS approach based on the large Fermi surface of Fermi liquid and band structure theory 
provides a good description of the overdoped region. At underdoping it is essential to include 
the emergence of a second energy scale, the pseudogap and its associated Gutzwiller factor, which 
accounts for a reduction in the coherent piece of the electronic Green's function due to increased 
correlations as the Mott insulating state is approached. In agreement with experiment, we find 
that the slope of the linear in T dependence of the low temperature specific heat rapidly increases 
above optimum doping while it is nearly constant below optimum. Our theoretical calculations also 
agree with recent data on Bi2Sr2_ x La x Cu06+,5 for which the normal state is accessed through the 
application of a large magnetic field. A quantum critical point is located at a doping slightly below 
optimum. 

PACS numbers: 74.72.-h,74.20.Mn,74..25.Bt 



I. INTRODUCTION 

The superconducting state of the underdoped cuprates 
shows anomalous properties when compared with the 
case of optimum or overdoped. There has been a con- 
siderable recent effort to understand these in a model 
which includes the emergence of a pseudogap for doping 
x below a quantum critical point (QCP) at x = x c . The 
model of Yang, Rice, and Zhang 1 (YRZ) on which this 
work is based, has the pseudogap forming about the anti- 
ferromagnetic Brillouin zone (AFBZ) boundary with its 
own characteristic energy scale. This model is different 
in many respects from other competing order proposals 
such as D-density waves 2 ' 3 - and from a preformed pair 
model 4 - which involves a single energy scale. For a review 
of the successes of the YRZ model in understanding the 
data, the reader is referred to Schachinger and Carbotte 5 . 
Many more details can be found in references^— It is 
important to understand that this work goes beyond ex- 
tensions of BCS theory to include effects such as strong 
coupling due to inelastic scatterin g 15-26 and possible 
anisotropy 2 ^— beyond a simple <i-wave superconducting 
gap. 

In a previous work^, we showed that the strong sup- 
pression of the specific heat jump at T c and correspond- 
ing reduction in condensation energy with increased un- 
derdoping can be understood as due to the emergence 
of a pseudogap. Here we consider the equally anoma- 
lous observation that the slope of the linear in temper- 
ature law as T — > is a strongly increasing function 
of doping x in the overdoped regime while it is nearly 
constant at underdoping^ 3 - - — While the overdoped case 
is characteristic of a Fermi liquid, the underdoped case 
requires a new ingredient for its understanding. In ad- 
dition, we will also consider the recent data of Wen et 



ali2£ on Bi2Sr2- x La x Cu06+,5 for which T c at optimum is 
only 30K and thus superconductivity can be quenched 
with a 9 Tesla magnetic field providing access to the nor- 
mal state as a function of doping in both Fermi liquid 
(overdoped) and pseudogap (underdoped) state. 

In Sec. II, we present the essential elements of the the- 
ory of YRZi needed for the calculation of the specific 
heat. Both pseudogap state alone and with the addition 
of superconductivity within a BCS formulation are con- 
sidered. For x greater than x c , the doping at which the 
QCP associated with the pseudogap formation is set, the 
Fermi surface is the usual large open contour of band 
theory and of a Fermi liquid. For x < x c , however, the 
Fermi surface reconstructs into Luttinger pockets which 
progressively shrink in size as the Mott insulating state 
is more closely approached. A consequence of this is 
that the density of states at the Fermi level is reduced. 
In addition, a Gutzwiller factor enters the theory which 
accounts for the depletion of the coherent part of the 
electronic Green's function due to increased correlations 
which eventually cause the transition to an insulating 
state. In Sec. Ill, we present our results for the spe- 
cific heat difference between superconducting and normal 
state as a function of temperature for values of doping 
ranging from 0.1 to 0.3 with x = x c = 0.2 the critical 
doping at optimum which is also the QCP in our model. 
The zero temperature limit of the specific heat difference 
is also considered more explicitly and compared with the 
density of states and the data of Wen et al^. In Sec. IV, 
we show that while pseudogap formation strongly affects 
the overall temperature dependence of the specific heat, 
it does not change its slope at low T. This arises be- 
cause this linear law only depends on the band structure 
near the Dirac point on the heavily weighted part of the 
Fermi surface in the nodal direction and this point is not 
changed by pseudogap formation in the model of Ref. 1 . 
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On the other hand the size of the slope depends directly 
on the Gutzwiller factor g t (x) which gives the magnitude 
of the remaining coherent part of the Green's function. 
It also depends on the ratio of the superconducting gap 
to critical temperature which is expected to strongly in- 
crease with decreasing doping x for x < x c .— These two 
effects combined lead to a rather constant value of the 
slope over a significant range of doping. For overdop- 
ing, Fermi liquid theory is recovered and in this case the 
slope shows much greater variation with x as in experi- 
ment. Comparison with data is presented. In Sec. V, we 
provide a summary and give our conclusions. 



II. FORMALISM 

In the resonating valence bond spin liquid model^, 
the self-energy due to the pseudogap is given by 
Ap g (fc)/[£°(fc)+w], where = — 2t(x)(cos k x a + cos k y a) 
with t(x) the effective first neighbor hopping parameter 
and Apg(fe) is a pseudogap which has ci-wave symmetry 
in the Brillouin zone. It is given by 

A (x) 

Apg(fc) = — — — (cosk x a — coskyd). (1) 

Here the amplitude Ap g (a;) is linear in doping x as shown 
in Fig. [IJa) where the superconducting dome is also 
shown for easy orientation. In the above, k is the mo- 
mentum and a is the lattice constant of the Cu02 plane. 
Without superconductivity, the coherent part of the elec- 
tronic Green's function is given by 



G(k,w) 



9t{x) 



c-e(fc)-A2 g (fc)/K°(fc) 



(2) 



where £(fc) is the electron dispersion curve of band the- 
ory and gt(x) is a Gutzwiller factor equal to 2x/ (1 + x). 
This latter quantity accounts for the effect of correlations 
which reduces the weight of the coherent part of G(k, uj) 
and adds an incoherent background not considered in 
this work. It also enters the renormalized band struc- 
ture dispersion curve which includes up to third nearest 
neighbor hopping and narrows as the Mott transition is 
approached with decreasing value of x. This narrowing is 
modeled by a second Gutzwiller factor g s (x) in addition 
to gt(x). Effectively, for a given k, there are two electron 
branches e£ with weights given by 



h k = o ± 



and 



wi = - 2 



1± 



[(a+^)/2] 2 + A2 g (fe) 



A2 g (fc) (3) 



(4) 



In terms of Ej~ , the Fermi surface contours of zero exci- 
tation energy are given by E^ = and these are shown 
in Fig. []Tb) for three values of doping x — 0.14, 0.18, and 
0.2. In the first two, there is a hole pocket centered about 
the nodal direction 6 — ir/4. This pocket is determined 
from the equation Ej~ = and Ej^ is positive only for mo- 
menta falling within the area defined by the hole pocket. 
The size of the pocket shrinks as x decreases and we come 
closer to half filling and the Mott insulating state. For the 
case x = 0.18, close to the QCP at x = 0.2 where pseudo- 
gap formation starts in our model, there is an additional 
electron pocket near the corner of the AFBZ. This pocket 
is determined by the equation E~£ = 0. Both electron and 
hole pockets have two sides, one weighted by Wj^ of order 
one and the other, which takes on a close resemblance to 
the AFBZ, has only a small weight in comparison. This 
small weight goes to zero in the limit of no pseudogap 
and the energy E^ becomes the Umklapp energy surface 
££. On the other hand, in this same limit, the heavily 
weighted part traces out the large Fermi surface of Fermi 
liquid theory (Fig.QJb), far right panel) which has weight 
one everywhere. These remarks make clear the evolution 
from large Fermi surface into small hole pockets. As the 
pockets shrink in size, the number of states which carry 
excitations of zero energy becomes small. Also, it should 
be kept in mind that because of the small weight on the 
backside of the hole pocket, we are effectively dealing 
with an arc when considering many properties. 

When superconductivity is included in a BCS formal- 
ism, the electronic spectral density A(k,u)) can be writ- 
ten in the form 



A(k,oj) = J2 9t{x)W^[{ulf5(u-El s ) + (vZ) 2 5(u- 



(5) 



with E% s = y/(E%) 2 + Aj c (k). Here, the supercon- 
ducting gap A sc (fc) is assumed to have the same d-wave 
dependence in momentum space as in Eq. (TTJ) for the 
pseudogap with amplitude A° c (x) replacing Ap g (a;) and 
Ag C (a;) is assumed to have the same doping dependence 
as the critical temperature T c (x) dome. This is shown in 
Fig.QJa) and for definiteness we will assume in all our nu- 
merical work that 2A® c (x)/kBT c (x) = 6. When later we 
compare with experimental data, we will remove this sim- 
plifying assumption. The Bogoliubov weights in Eq. ([5]) 
are 
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id 



E« 



E a s 



(6) 
(7) 



The dispersion curves for and are taken from the 
work of Refi and are unchanged here as are the other 
parameters, namely, 



A° c (» = 0.14to [1 - 82.6(a; - 0.2) 2 ] 



(8) 
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FIG. 1: (Color online) (a) Phase diagram for the model dis- 
cussed in the paper, showing the pseudogap Ap g and super- 
conducting gap A° c in units of to as a function of doping x. 
(b) Illustration of the reconstruction of the Fermi surface in 
the first quadrant of the Brillouin zone for doping values of 
x = 0.14, 0.18 and 0.2. The red dashed line is the AFBZ 
boundary. 



and 



(9) 



with to an unrenormalized nearest neighbor hopping pa- 
rameter characteristic of the Cu02 plane. 

There are two equivalent ways for calculating the spe- 
cific heat Cy{T) = ^{T)T. One is through the entropy 
S(T), the other through the internal energy U(T). The 
entropy is 

S(T) = -2k B g t (x) J2 WZ{f(E% s )lnf(E% s ) 

k,a— ± 

+[!-/»)] Mi -/(*&)]}, 

(10) 

where f(x) is the Fermi-Dirac temperature distribution 
function. Alternatively, the internal energy can be ex- 
pressed in terms of the single spin density of states N(cu) 
as 



U(T) = 2 / ujN{uj)f(uj)duj 



(11) 



where the two is for spin degeneracy and the density of 
states is given by 

N(cj) = J2 9t(x)Wn(uZ)H(oj-El s )+(v%) 2 6(u+El s )}. 

(12) 

In both Eqs. (flU)) and (fl"2"l) , the sum is over the entire 
Brillouin zone and C V (T) = dU(T)/dT = j(T)T. 
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FIG. 2: (Color online) Aj(T) versus T/t , where A7 = -y sc - 
j„ is in units of k%/to per volume. Curves are shown for 
(a) underdoped and (b) overdoped cases along with optimal 
doping for reference. 



III. NUMERICAL RESULTS FOR THE 
SPECIFIC HEAT 

In Fig. [21 we present our results for the difference be- 
tween the superconducting and normal state specific heat 
7(T), namely A'y(x, T) = 7 sc (a:, T) - J n (x, T) in units of 
k\ /to per volume as a function of temperature T. Several 
values of doping x are considered as indicated in the fig- 
ure. The top frame covers the underdoped regime, while 
the bottom frame is for overdoped. As the superconduct- 
ing dome given in Eq. is symmetric in doping about 
optimum x = 0.2, the top and bottom frame curves come 
in pairs with the same value of critical temperature (with 
the exception of not displaying x = 0.3 for clarity). Com- 
parison of top and bottom curves in a given pair shows 
that the formation of a pseudogap and associated Fermi 
surface reconstruction [Fig. QJb)] provides a drastic sup- 
pression of the jump at the critical temperature T c and 
also of the slope just below T c . This is in qualitative 
agreement with experiment— ~—, as discussed already in 
RefA These results are clearly not part of ordinary d- 
wave BCS theory where instead the jump is large and 
relatively independent of doping as seen in the Fig.[^b). 
For reference in assessing these curves, we recall that for 
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FIG. 3: (Color online) (a) A7 as T — > normalized to the 
same quantity at optimal doping versus doping x. Values 
extracted from the data of Wen et al.— are shown along with 
the DOS at zero energy N(0) similarly normalized. The inset 
shows N(uj) in units of to for several dopings as a function 
of w/to- N(lj) and N(0) refer to the normal state with a 
pseudogap. (b) Value of A7 at T c from Fig. [2] (left axis) and 
slope just below T c , d/l\'y{T c ) / dT (right axis) as a function of 
x. 



a constant density of states model with a superconduct- 
ing d-wave gap defined on the Fermi surface, the canon- 
ical value of the jump A-j(T c )/^ n (T c ) is 0.95 for a gap 
to critical ratio 2Ag C /fcsT c = 4.3 and is the same for all 
superconductors regardless of the size of T c . Here it de- 
viates from this universal law because our superconduct- 
ing gap is defined in the entire Brillouin zone according 
to Eq. ((T|) with the pseudogap replaced by A° c and our 
energy bands can be complicated even when the large 
Fermi surface of Fermi liquid theory is involved. Also, 
the magnitude of the jump itself is increased because we 
have used a gap to T c ratio of 6 rather than the weak 
coupling limit of 4.3 and this has resulted in a change 
of A 7 (T C )/7„(T C ) from 0.95 to - 1.7. Thus, the large 
decrease in the jump A7(T C ) seen in Fig. HJa) is a direct 
consequence of pseudogap formation and accompanying 
Fermi surface reconstruction. 

The size of the specific heat difference seen at zero 
temperature is also of interest and is very different in 



the underdoped than in the overdoped regime. While for 
near optimum, optimum and overdoped cases the value of 
A 7 (T) as T — > is around -1 in our units, for the under- 
doped cases, it has moved instead to a value of roughly 
~ —0.5. In this limit Aj(T — 0) simply reduces to its 
normal state value 7™ and is a direct measure of the value 
of the density of states N(oj) at the Fermi energy, i.e. at 
uj = 0. In Fig.[3]Ja), we compare the results of our calcu- 
lations for A7(T = 0)/A7opt(7 1 = 0) (circles) with the 
specific heat data of Wen et ali^S on Bi 2 Sr2_ x La x Cu06+5 
(squares). Here, OPT refers to the value taken at opti- 
mal doping (x = 0.2 in our model). In carrying out this 
comparison, we have fit a parabolic form to the data of 
Wen et al. for T c versus doping. We found that the fit 
which best captured the data was when optimum was 
taken at their doping of p = 0.165 and to is set to 56 
meV for 2Ao//csT c = 6. Similar, but less ideal fits would 
also work for optimum at p — 0.16. Note that this ma- 
terial does not have a very high T c and consequently has 
a narrower dome compared with what is typically found 
in the literature. It was chosen for experiment because 
reasonable magnetic fields can be used to suppress T c to 
zero and so access the normal state. Given our fit, we 
have shifted the value of p in experiment by 0.035 to get 
correspondence with our theoretical work for which x at 
optimum is kept at x = 0.2. Further discussion of this is 
presented in Sec. IV. We have kept the pseudogap line to 
be as shown in Fig. [TJa) which coincides with optimum 
doping, but this could easily be changed. Returning to 
Fig. EJa), both theory and experiment show two distinct 
regimes: a somewhat flat or slowly rising region above 
x ~ 0.17 with a sharp drop below this value. The behav- 
ior in the region above x = 0.2 is sensitive to the choice of 
bandstructure parameters as shown in the original work 
of YRZ (Rcf. 1 ). Because the data is fairly flat in this 
region, we have opted to alter the bandstructure accord- 
ingly and have allowed the t' and t" parameters to con- 
tinue to vary with x rather than to become constant for 
x > 0.2 as in Refi Returning to the rapid drop around 
x ~ 0.17, we find that it is displaced to slightly lower 
values of x in the data as compared to theory and we 
take this to mean that in Bi 2 Sr2_ x La x Cu06+5 the QCP 
associated with the start of pseudogap formation occurs 
at a doping level close but slightly below optimum. It is 
important to realize that in the present model, the point 
at which the value of the density of states at the Fermi 
level N(0) begins to be strongly reduced by the growth 
of the pseudogap, is not exactly at the QCP but is in- 
stead displaced to smaller x values. The physics of this 
displacement is easily understood with reference to the 
inset in Fig. [3fa). First we show, in the main frame, 
that our numerical results for N(0) versus x (triangles) 
are very close to the solid black circles as they must be. 
These were obtained not from an extrapolation to zero 
temperature of a specific heat calculation as were the 
solid circles but instead from the full DOS N(u>) versus 
uj shown in the inset. This good agreement is taken as 
providing a check on our numerical work. More signif- 
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icantly it is important to realize that the Fermi surface 
reconstruction from the large open Fermi liquid surface 
of Fig. [T] (x = 0.2) to a Luttinger hole pocket does not 
immediately lead to change in the value of N(0). In fact, 
as can be seen clearly in the inset of Fig.[3ta), at x = 0.18 
the DOS at u = is hardly changed from its Fermi liquid 
value (see solid black curve for x = 0.2 for comparison). 
Rather the effect of the pseudogap is to provide a de- 
pletion of states at negative w. For x = 0.16 however, 
the upper edge of this depletion region which is rather 
sharp, has moved across the Fermi energy and N(0) has 
become significantly reduced as shown in the triangles 
of the main frame in Fig. [3][a). This is not surprising. 
The DOS at cj — depends on the Luttinger contours 
of zero energy shown in Fig. [T] At x = 0.18, there are 
both hole and electron pockets and the number of states 
having zero energy is not very different from the num- 
ber when the large open Fermi surface applies. It is only 
when we reach a single small hole pocket and no electron 
pockets that depletion of states at cj = becomes signif- 
icant. On the face of it, one might even think that more 
zero energy states are involved when x — 0.18 than when 
x = 0.2 because the length of the Fermi surface is larger 
when there are pockets. But this is not so because, as 
we have already stated, parts of the Luttinger contours 
carry very little weight. 

In Fig. [3][b), we show results for the jump at T c of 
the specific heat A7(T C ) (circles) and its slope out of T c 
(squares) as a function of doping. We note that for the 
jump AjfTc), a notable drop occurs almost immediately 
below optimum x — 0.2, which can be used to identify the 
QCP. For the slope, the signature of the QCP is not as 
sharp, only a change in curvature arises at x — 0.2. These 
results are to be contrasted with those in Fig. EJa). If 
one were to estimate the value of doping to be associated 
with the QCP from the behavior of A 7 (T) at T -> 0, 
it is necessary to account for a significant displacement 
downward towards smaller values of x of the rapid drop 
in N(Q) as compared to the value of x at the QCP. Thus 
extracting a QCP from thermodynamics requires some 
care but it clearly can be done either from the jump at 
T c or the value of the difference Aj(T) as T — >• 0. 



cone is the only region in momentum space where there 
is a finite thermal occupation of excited quasiparticles. 
Thus, the specific heat can depend only on characteris- 
tic parameters associated with the Dirac point. But this 
point is particularly simple. It corresponds to ET = 
in the normal pseudogap case and does not change with 
the onset of superconductivity. In the nodal direction, 
the pseudogap A pg (fc) = and therefore = which 
is the condition for the underlying large Fermi surface. 
Note that strictly speaking involves the chemical po- 
tential fi p associated with the case when the pseudogap 
is present and is slightly different from the chemical po- 
tential of Fermi liquid theory. In the former case, it is 
determined from the Luttinger sum rule while in the lat- 
ter case of a large Fermi surface one might determine 
\x v from the DOS filling up to k = kp. Neglecting this 
small difference, the band energy at the Dirac point is un- 
changed from its Fermi liquid value and the well known 
techniques 37 for obtaining the u> — > limit of the DOS 
N(uj) apply unaltered. The result is 



N(w) 



I 9t(x) 
nvpVA 



(13) 



where vf and va are the Fermi and gap velocity, respec- 
tively. Note that Eq. ([13")) predicts that the only effect of 
the pseudogap formation on the DOS around the Dirac 
point is the appearance of the Gutzwiller factor gt(x). 
While to a good approximation vp is unchanged, the gap 
velocity can be changed in magnitude if the ratio of zero 
temperature gap to critical temperature T c is affected by 
pseudogap formation as was found in the recent work of 
Schachinger and Carbotte^. This provides a second im- 
portant change in Eq. (| 13[) as compared with the more 
familiar Fermi liquid case. 

Inserting Eq. (fT3"|) into Eq. (|11[) for the internal energy 
U(T), we obtain the simple analytic result for 7 SC (T) in 
the limit of low temperature: 



lsc (T) = Ak%(k B T)-^-h, 



irvpVA 



(14) 



where k B is the Boltzmann constant and h is a number 
given by 



IV. LOW TEMPERATURE BEHAVIOR 

We now consider the low temperature behavior of the 
superconducting state specific heat 7 sc (r) which is em- 
phasized in Fig. Eta)-(d) for underdoped cases. We base 
our discussion on the physics of the schematic shown in 
Fig. IDJe). The Luttinger pockets are shown as heavy 
black lines for x = 0.18 which has electron as well as hole 
areas. Also shown is a Dirac cone centered on the highly 
weighted Fermi surface contour in the nodal direction. 
The cone illustrates the quasiparticle energies at low en- 
ergy in the superconducting state as a function of k x and 
k y in the Brillouin zone. Only the upper right quadrant 
is depicted. At very low temperatures, the tip of this 



h = A 



dy|?/| 3 cosh 2 (y) ~ 5.4. 



(15) 



We have obtained, as in ordinary d-wave BCS theory, a 
linear in T law with the same material factors appear- 
ing, i.e. vf and va, but with an extra Gutzwiller fac- 
tor of gt(x). We also need to note that A can be af- 
fected by the variation with doping x of the gap ratio 
R(x) = 2A° sc (x)/k B T c (x). 

In Fig. |H we show our numerical results for 7 SC (T 1 ) at 
four values of x namely 0.12, 0.14, 0.16 and 0.18. In all 
cases, the solid line is for the case when no pseudogap is 
included while the dashed-dotted includes a pseudogap 
with its magnitude chosen to correspond to the value of 
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FIG. 4: (Color online) j(T) versus T '/to shown in the absence 
of any gaps jm, with superconductivity only jsc, pseudo- 
gap only jpg and with both gaps present jsc+pg- The red 
dashed line gives the theoretical expression for the slope from 
Eq. (|14[1 . (a)-(d) show a range of dopings for the underdoped 
case, (e) Schematic of the superconducting Dirac cone shown 
on the side of the x = 0.18 Fermi pocket with large quasipar- 
ticle weight. 



doping x chosen. Otherwise, there are no other changes. 
Gutzwiller factors are included in these curves and the ra- 
tio R(x) is fixed at value 6 as in all other numerical work 
presented in this paper. While the inclusion of a pseu- 
dogap has drastic effects on j sc (T), the slope as T — > 
is completely unaffected and this slope agrees perfectly 
with the simplified analytic results of Eq. ([14")) given as 
the red dashed line. This constitutes an important pre- 
diction of YRZ theory and will be verified later when we 
make comparison with experimental data. 

The heavy dotted lines in Fig. 2] show the normal state 
specific heat ^^(T) in the Fermi liquid and the light dot- 
ted, with a pseudogap included jn+pg{T). It is instruc- 
tive to consider these in some detail and in particular to 
describe how they are related to the detailed variation 
of N(uj) versus lu in the pure pseudogap state, i.e. with 
no superconductivity. We begin with the case x = 0.18 



for which the pseudogap is small and the Fermi contours 
include electron as well as hole pockets. In this case, 
the light dotted curve falls below the heavy dotted one 
but rises to meet it as T — > 0. This behavior can be 
traced to the lu variation seen in the inset of Fig. [3ta) 
which shows N(lu) versus lu. The solid black curve forms 
a reference and is the Fermi liquid result for x = 0.2, 
zero pseudogap. Comparing with the long-dashed green 
curve for x — 0.18, we note no visible change in the value 
of N(lu — 0) but there is a significant depression of DOS 
at negative energies. As T — ► only N(Q) is sampled in 
the specific heat and hence both cases, with and with- 
out a pseudogap, agree. At small but finite T, however, 
the light-dotted curve falls below the heavy dotted line 
in Fig. 0Jd) because the dip in its DOS at lu < starts 
to be sampled and this reduces the specific heat. Specific 
heat is, however, a rather broad spectroscopy for N(lo) 
because the thermal factor in j(T) in Eq. (fTTj) samples 
of order 5&bT or so about lu = 0. The case x = 0.12 
is also noteworthy. In this instance the DOS about the 
Fermi energy is nearly monotonic although depressed in 
value as compared to x = 0.2 by about 50% and this 
correspondingly reduces the value of jn(T) by the same 
amount with little other changes. 

We turn next to a comparison with data on the slope of 
the specific heat in the T — !> limit. These are presented 
in Fig. [5fa) for the Bi2Sr2- x La x Cu06+<5 data of Wen et 
al^ and in (b) for the data on Bi-2212 and (Y,Ca)-123 of 
Loram et ali^r— where we have read slopes from the pub- 
lished figures as best as was possible. A more thorough 
analysis by the experimentalists might achieve better ac- 
curacy. In the inset of Fig. (5]Ja), we show our results for 
the fit to the T c dome of the Wen et ali^ data. We find 



225(p- 0.165) 2 ]i , 



(16) 



which provides a good representation of the T c data as 
a function of doping p. We then take x = p + 0.035 to 
place the data on our curves. In the main frame, the 
solid black circles are the data for the slope at T — » 
normalized to its value at optimum doping. The solid 
red line gives results of Eq. (fT^Ij) with the assumption 
of R = 6. The agreement is very good and the sharp 
rise in slope in the overdoped region is captured by our 
model. For the deeply underdoped case, our theoreti- 
cal values are somewhat higher than experiment, but as 
we have already mentioned, the gap ratio may well vary 
with doping. Solving a BCS-like pairing equation with 
pseudogap formation and corresponding Fermi surface re- 
construction accounted for, Schachinger and Carbotte 5 
found that (for x < 0.2) approximately 



R(x) -4.3[l + 75(a;-0.2) 2 ]. 



(17) 



If this correction, scaled to 6 at optimum, is incorpo- 
rated into the comparison with experiment, we get the 
red dashed curve which agrees better with the data at 
small x. Note finally that to compare data with theory, 
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FIG. 5: (Color online) The slope of the specific heat at low T 
normalized to the value at optimal doping. Frame (a) shows 
the data extracted from Wen et al.— compared with the the- 
ory of Eq. (|14|) for the T c dome given in the inset where we 
have fit a functional form for x to their T c data (see text for 
discussion). The dashed curve includes the correction R(x) 
discussed in the text, (b) The extracted data of Loram et 
al.— ~_3£ j s compared to the theoretical slope. Shown also is 
the slope data from the STS DOS of Pushp et al.— . 



we have shifted all x values in Fig.[5ja) by 0.035 because 
we wished to remain, as in the paper of YRZ 1 , with op- 
timum doping at x = 0.2 rather than the experimen- 
tal value of about 0.165. This is also true for Fig. [SJb) 
where we compare with data on Bi2212 (triangles) and 
(Y,Ca)-123 (squares). The agreement with theory (solid 
and dashed red curves) is again good. We have also in- 
cluded one further comparison with scanning tunneling 
spectroscopy (STS) results on Bi2212 by Pushp et alM 
(circles). STS gives the DOS and not the specific heat, 
but this latter quantity follows directly from a knowledge 
of N(ia). We can use the STS data in the limit of ui — >• 
to determine N(u) as in Eq. (H"3l) and so get j S c(T) for 
T — >• 0. These are the solid circles which provide a con- 
firmation of the specific heat results and also provide a 
significant cross check between these two important but 



very different probes of the microscopic structure of the 
superconducting state in the underdoped cuprates and 
indeed over the entire phase diagram. 



V. SUMMARY AND CONCLUSIONS 

The behavior of the specific heat of the underdoped 
cuprates differs profoundly from that observed on the 
overdoped side of their phase diagram. At optimum 
and overdoping, an ordinary BCS approach based on 
a Fermi liquid normal state with constant DOS pro- 
vides a first reasonable understanding on the assump- 
tion that the gap has ci-wave symmetry. This ensures 
a linear in T low temperature law for the specific heat 
7 (T) with CV(T) = 7(T)T. It also gives a jump at 
T c normalized to its normal state value of 0.95 and this 
value can be increased if the gap to T c ratio is changed 
from 4.3 to a higher value. On the deeply underdoped 
side, however, there is new physics which cannot be de- 
scribed even when strong coupling effects^— due to in- 
elastic scatteringi2r— are accounted for and/or effects of 
anisotropy2£~— . which go beyond a simple lowest har- 
monic e?-wave picture for the superconducting gap. 

In the resonating valence bond spin liquid picture, a 
second energy scale, the pseudogap, emerges and grows 
in magnitude as x is reduced towards the Mott insulating 
state. This pseudogap leads to a loss of metallicity. It 
also radically reduces the size of the specific heat jump 
at T c as noted in experiment^— and also in theory^. 
However, as we show here the limit of T — > of 7 SC (? 1 ) 
is not directly affected by the size of A° g because this 
limiting value depends only on the band structure and 
superconducting gap right at the Dirac point in the Bril- 
louin zone. But this point is not importantly changed 
by pseudogap formation and Fermi surface reconstruc- 
tion so that the formula for the slope dj sc (T)/dT as 
T — > remains unchanged in form but with two im- 
portant modifications. First, in the resonating valence 
bond spin liquid, there appears a Gutzwiller factor gt[x) 
which depends strongly on doping x and represents the 
remaining weight in the coherent part of the electronic 
Green's function as correlations become more important 
and consequently shift more spectral weight into an in- 
coherent background at higher energies. A second factor 
is that the ratio of 2A° c (x)/kBT c (x) = R(x) can vary 
with doping. For overdoped and optimally doped, in or- 
dinary BCS d-wave it has a value of 4.3, but Schachinger 
and Carbotte^ have found that it rises considerably in 
the underdoped region of the phase diagram. These au- 
thors solve a BCS gap equation generalized to include 
pseudogap formation and attendant Fermi surface recon- 
struction. The results of such a theoretical study show 
that R(x) increases beyond a value of 7 before the lower 
end of the superconducting dome is reached. These two 
effects, along with the linear in T law which we have 
shown to still hold in YRZ theory at underdoping, al- 
low us to understand a previously anomalous feature of 
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the data namely the slope remains reasonably constant 
at undcrdoping while it increases sharply in the over- 
doped region. This last observation is consistent with 
a slope which varies inversely as the gap and this gap 
decreases toward zero as we approach the upper end of 
the dome. Comparison between theory and data on BiS- 
rLaCuO, Bi2212, and (Y,Ca)-123 show good agreement. 
A further comparison was made with STS data which 
provides information on the average density of quasipar- 
ticle states. It was noticed by Pushp et alM- that the 
slope of this quantity, in the lo — > limit, while increas- 
ing with reduced value of T c at overdoping, saturates and 
perhaps even decreases slightly with decreasing x in the 
highly underdoped regime. But the low lo dependence of 
the DOS determines the low temperature behavior of the 
specific heat. For the specific case of Bi2212, we found 
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good agreement between STS and specific heat data fur- 
ther confirming our work and providing a strong test of 
the consistency between experimental data obtained by 
these two very different techniques and their consistency 
with the resonating valence bond spin liquid. 
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